Low-Density Aerodynamics of the Stardust Sample Return Capsule 
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The aerodynamics of the Stardust Sample Return Capsule are analyzed in the low- 
density, transitional flow regime using free-molecular, Direct Simulation Monte Carlo, 
Navier-Stokes, and Newtonian methods to provide inputs for constructing a transitional 
flow bridging relation. The accuracy of this bridging relation in reconstructing the aero- 
dynamic coefficients given by the more exact methods is presented for a range of flight 
conditions and vehicle attitudes. There is good agreement between the various prediction 
methods, and a simple sine-squared bridging relation is shown to provide a reasonably 
good description of the axial force, normal force, and pitching moment over a range of 
Knudsen numbers from 0.001 to 10. The predictions show a static instability of the Star- 
dust capsule in the free-molecular regime that persists well into the transitional flow. 
The addition of a thin disk to the base of the capsule is shown to remove this static 
instability. However, the extremely high entry velocity of 12.6 km/s for the proposed 
trajectory introduces difficult design issues for incorporating this disk caused by the high 
aerothermal loads that occur even under relatively rarefied conditions. 


Nomenclature 

A re j Reference area (= ttRC), m 2 

c.g. Center of gravity (= 0. 36751b) , m 
Ca Axial force coefficient, Fa/ (\pooVoo 2 A re j) 

Cn Normal force coefficient, Fn/ {\pcaVcC 2 A re j) 

C' m Pitching moment coefficient, 

My/{\ 

poo too 2 A re j L re j) 

Fa Axial force, N 

Fn Normal force, N 

Db Maximum body diameter, m 

Kn Knudsen number, A / Db 

L re f Reference length (= Db), m 

M Mach number 

M y Moment about pitch axis, N-m 

n Number density, 1/m 3 

P Pressure, N/m 2 

q Heat flux, W/m 2 

R a Afterbody base radius, m 

Rb Maximum body radius, m 

R n Nose radius, m 

R s Shoulder radius, m 

S Distance along body, m 

T Temperature, K 

V Velocity, m/s 

x, y, z Cartesian body axes, m 

a Angle of attack, deg 

A Mean free path, m 
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p Mass density, kg/m 3 

r Shear stress, N/m 2 

9 a Afterbody half-cone angle, deg 

Of Forebody half-cone angle, deg 

Subscripts 

eg Center of gravity 

cont Continuum 

fm Free molecular 

oo Freestream 

Introduction 

S TARDUST is a NASA Discovery-class mission de- 
signed to make a close encounter with the comet 
Wild-2, collect samples of dust and volatiles within the 
coma and return them to Earth. 1 The Earth return 
will use a trajectory that produces an entry velocity 
of approximately 12.6 km/s. Landing footprint restric- 
tions require accurate knowledge of the aerodynamics 
of the reentry capsule throughout the trajectory from 
the free-molecular flow regime of initial atmospheric 
encounter to the subsonic regime just prior to landing. 
The high entry velocity makes the size of the landing 
footprint somewhat more sensitive to the low-density 
aerodynamics than typical reentries from low Earth 
orbit. More importantly, aerodynamic instabilities in 
the low-density regime can cause angle-of-attack ex- 
cursions that produce undesirable capsule attitudes at 
peak heating conditions. 

The low-density aerodynamics of entry vehicles 
are frequently determined through the use of free- 
molecular and Newtonian methods with bridging rela- 
tions to define the transitional flow aerodynamics. 2,3 
These bridging relations typically use the Knudsen 
number as the independent parameter and require 
knowledge of the Knudsen number above which molec- 
ular collisions can be neglected (free-molecular) and 
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the Knudsen number below which the flow may be 
treated as a continuum. The aerodynamic coefficients 
in the free-molecular and continuum limits are also 
required. In many instances, this approach produces 
estimates of the aerodynamic coefficients that are suffi- 
ciently accurate especially when combined with active 
or passive flight control that has sufficient margins 
for correcting entry trajectory and attitude. How- 
ever, because of the need for increased accuracy for the 
passively-controlled Stardust Sample Return Capsule 
(SRC), a study was conducted in which both direct 
simulation Monte Carlo (DSMC) and Navier-Stokes 
calculations were used to accurately define the transi- 
tional aerodynamics. 

The purpose of this paper is to discuss the method- 
ology and present the results used to define the low- 
density aerodynamic database for the SRC. A descrip- 
tion of the methodology used to define the complete 
aerodynamic database from hypersonic, free-molecular 
to subsonic, continuum flow is given in a companion 
paper. 4 Aerodynamic drag based on free-molecular, 
2D and 3D DSMC, Newtonian, and 2D and 3D Navier- 
Stokes methods is used to establish the Knudsen num- 
ber limits of the transitional flow regime and thereby 
the appropriate bridging relation parameters. Further 
3D DSMC analyses are then used to establish the va- 
lidity of the bridging relations for predicting other 
aerodynamic coefficients including the pitching mo- 
ments needed to determine static stability. Finally, 
since the baseline configuration exhibits a static in- 
stability, additional free-molecular and DSMC results 
are shown for a modified geometry representing one 
possible remedy to the static instability. 


Geometry 

The baseline geometry of the Stardust Sample Re- 
turn Capsule is shown in Figure 1 and consists of a 
60-deg half-angle, spherically-blunted cone forebody 
with a 30-deg half-angle conical afterbody. The nose 
radius, R n , is 0.2286 nr, the shoulder radius, R s , is 
0.01905 nr, and the overall diameter, Db, is 0.8128 nr. 
The afterbody base has a radius of 0.2116 nr. The cen- 
ter of gravity of the baseline configuration is 0.2987 nr 
aft of the nose ( x cg /Db = 0.3675). This aft location 
results in a static instability in the free-nrolecular flow 
regime. One approach to removing this instability 
without moving the center of gravity is analyzed which 
has a thin circular disk or “skirt” added to the after- 
body base as shown in Figure 2. The diameter of this 
disk is the same as the overall diameter of the fore- 
body, 0.8128 nr. The thickness of the disk was not set 
by the design at the time of analysis, but for the pur- 
pose of this analysis was assumed to have a thickness 
of 0.0015 nr. 



Fig. 1 Baseline geometry of Stardust Sample Re- 
turn Capsule. 


Skirt 



Fig. 2 Modified geometry with afterbody skirt. 

Trajectory/Flight Conditions 

A nominal reentry trajectory was determined from 
3 degree-of-freedonr analyses, 4 and specific trajectory 
points selected for detailed analysis are given in Ta- 
ble 1. Because of the small size of the capsule, atmo- 
spheric interface denoting the beginning of the transi- 
tional flow regime occurs at a relatively low altitude. 
For this trajectory, free-molecule flow is expected at 
altitudes above about. 130 km, and fully continuum 
flow is expected at altitudes below about. 80 km. Ax- 
isynnnetric DSMC calculations were performed for all 
trajectory points given in Table 1 at zero degrees angle 
of attack. Full 3D DSMC calculations were performed 
for all points except, the lowest, altitude at zero de- 
grees and for altitudes of 134.75, 100.9 and 83.68 km 
at angles of attack of 10 and 30 deg. Axisynmret.- 
ric Navier-Stokes calculations were performed for an 
altitude of 83.68 km at 0 deg. and full 3D Navier- 
Stokes calculations were performed for angles of attack 
of 0 and 10 deg. Additional Navier-Stokes calculations 
were performed at several altitudes lower than those 
in the table for development, of the complete aerody- 
namic database. 4 
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Table 1 Flight Conditions for Entry Simulations 


Altitude, 

Velocity, 

Number 

Mole Fractions 

T 

- L 00 5 

T 

- L W 5 

Knudsen 

Mach 

km 

m/s 

Density, 1/m 3 * 

0 2 

n 2 

O 

K 

K 

Number 

Number 

134.75 

12,597.1 

1.48832xl0 17 

0.0659 

0.6716 

0.2625 

577.23 

400 

12.8 

23.96 

120.45 

12,607.7 

5.93941xl0 17 

0.0845 

0.7327 

0.1828 

381.15 

500 

2.92 

30.27 

100.90 

12,620.2 

1.09988xl0 19 

0.1768 

0.7844 

0.0388 

199.37 

1,000 

0.136 

43.94 

92.00 

12,618.5 

4.98474xl0 19 

0.2056 

0.7873 

0.0071 

202.05 

1,200 

0.0301 

44.10 

83.68 

12,591.5 

1.77978xl0 20 

0.2385 

0.7615 

0.0000 

216.54 

1,500 

0.00857 

42.68 

75.98 

12,486.8 

5.73854xl0 20 

0.2385 

0.7615 

0.0000 

218.14 

1,800 

0.00266 

42.17 


Computational Methods 

Bridging Relation 

The form of the bridging relation used in this study 
is: 

C'A — CA,cont T {CaJui CJA,cont ) SZU (</) (1) 

where 

<f> = 7T (cq + a 2 log w Kn) (2) 

Here, C'a represents the axial force coefficient, but 
similar expressions are used for other aerodynamic co- 
efficients. The quantities C'A,cont and Cajhi represent 
the values of the coefficient in the continuum and free- 
molecular limits respectively. The constants a i and 02 
are determined by choosing Knudsen numbers corre- 
sponding to each of these limits. While general guide- 
lines frequently give these limits to be Kn=0.01 for the 
continuum limit and Kn=10.0 for the free-molecular 
limit, it has been shown that both of these limits fre- 
quently need to be extended to achieve fully continuum 
or fully free-molecular flow. 2,5 Furthermore, since the 
constants, a 1 and a 2 , are simply adjustable parame- 
ters in Eqn. 1, they should be adjusted to give the 
best overall description of the transitional flow aero- 
dynamics when additional data are available. 

In this study, free-molecular, DSMC, Navier-Stokes, 
and Newtonian calculations are used to define the 
behavior of the axial-force coefficient throughout the 
transitional regime. The constants, a 1 and a 2 , are 
determined to give the best overall match to the com- 
puted C'a at a = 0°. Continuum and free-molecular 
values are then computed for all aerodynamic quan- 
tities (axial force, normal force, and pitching moment 
coefficients) over a range of a and used together with 
the constants a 1 and a 2 to provide a complete aero- 
dynamic database. To verify the adequacy of the 
bridging relations for non-zero incidence angles, full 
3D DSMC calculations are performed for the selected 
incidence angles and trajectory points described pre- 
viously. 

Free-Molecular and DSMC Analyses 

Free-molecular analyses were initially performed 
for the baseline configuration using standard free- 
molecular aerodynamic relations applied to a dis- 
cretized surface geometry. However, the concave shape 


of the modified geometry that later needed to be ad- 
dressed made the use of typical line-of-sight shadowing 
somewhat suspect. Therefore, the results denoted free- 
molecular in this paper were actually calculated using 
DSMC with molecular collisions disabled. Compar- 
isons of the free-molecular results with the collisionless 
DSMC for the baseline configuration (which has no 
concavities) show that the two methods give aerody- 
namic coefficients that agree to within less than 1%. 

DSMC analyses used the standard algorithm devel- 
oped by Bird. 6 Axisymmetric analyses were performed 
using the standard G2 code 7 and full 3D analyses 
were performed using a relatively new code called 
DAC developed by LeBeau. 8 The 3D code uses an 
unstructured triangular grid to represent the surface 
geometry together with a two-level embedded Carte- 
sian grid for the flowheld. Gas collisions are modeled 
using the variable hard-sphere (VHS) model and the 
Larsen-Borgnakke model for internal energy exchange. 
Gas-surface interactions are assumed to be fully dif- 
fuse with full energy accommodation to the surface 
temperature. The surface temperature is based on 
an approximate heat transfer estimate and a simple 
radiative-equilibrium analysis and is given in Table 1 
for each of the trajectory points analyzed. In gen- 
eral, both the axisymmetric and 3D analyses model 
the gas as 5-species reacting air using a 23-equation 
reaction set. However, for the lowest altitude cases 
with the 3D code, it was necessary to use a parallel 
version of the DAC code in which chemistry is not yet 
implemented. Therefore, some limited comparisons 
were made between 3D calculations with and without 
chemistry based on relatively coarse grids for angles 
of attack up to 30 deg. These comparisons show that 
while there is considerable dissociation of O 2 and N 2 
for both the 100.9 and 83.68 km cases, this dissociation 
has a negligible effect on the aerodynamics. Therefore, 
the 3D results presented at 83.68 km were calculated 
without considering chemical reactions. 

Navier-Stokes and Newtonian Analyses 

Navier-Stokes analyses are performed using the Lan- 
gley Aerothermodynamics Upwind Relaxation Algo- 
rithm (LAURA) CFD code . 9,10 LAURA is an upwind- 
bias, point-implicit relaxation algorithm for obtaining 
the numerical solution to the Navier-Stokes equations 
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for three dimensional, viscous, hypersonic flows in 
thermochemical nonequilibrium. It has been used to 
describe the aerodynamics of several blunt bodies sim- 
ilar to the Stardust Sample Return Capsule including 
Mars Pathfinder. 11 Calculations were performed at the 
trajectory point corresponding to the lowest altitude 
given in Table 1 (as well as several lower altitudes 
given in Ref. 4 for the continuum portion of the trajec- 
tory) with no-slip boundary conditions imposed. The 
LAURA results are used to establish the magnitude 
of the aerodynamic coefficients at. the Knudsen num- 
ber corresponding to the continuum limit. Newtonian 
analysis is then used to provide the variations of these 
coefficients with angle of attack relative to the LAURA 
results. 

Results 

Flowfields 

Number density contours for three different trajec- 
tory points at a = 30° are shown in Figure 3 to 
illustrate the changing character of the flowfield as the 
degree of rarefaction changes. At 135 km, the flow is 
essentially free molecular and while there is a signifi- 
cant. increase in density in front, of the vehicle, there is 
no distinct, shock behavior. As the density increases 
at. the lower altitudes, a. shock layer becomes more 
evident, until at. 83 km, a. very distinct, bow shock is 
produced, and the flow does not. expand in the wake 
t.o densities as low relative t.o the freest.rea.m as it. does 
at. the higher altitudes. 

Baseline Geometry Aerodynamics 

Figure 4 shows a. comparison of the axial force coef- 
ficient. computed by the various methods as a. function 
of Knudsen number for a = 0°. There is generally 
good agreement, between the 2D (a.xisymmet.ric) and 
3D DSMC methods as well as the Na.vier-St.okes results 
at Kn=0. 00857. The bridging relation is found to give 
a. reasonable fit. t.o the computed results by choosing 
Kn=0.001 as the continuum limit, and Kn=10 as the 
free-molecular limit.. The lower limit, is reduced by 
a. factor of 10 below that, value normally considered 
to represent, the start, of continuum flow. The DSMC 
results indicate that, the free-molecular limit, should 
perhaps be increased, but. it. was felt. that, more weight, 
should be given t.o the lower Knudsen numbers where 
the aerodynamic forces are greater and have more in- 
fluence on the trajectory. Other forms of bridging 
relations might, also give a. better fit. t.o the data.. 2,5 

One reason the predicted results do not. follow the 
current, bridging relation better is that, the points cho- 
sen follow a. computed trajectory where the vehicle 
velocity and atmospheric properties are changing as 
well as the atmospheric density. Furthermore, the sur- 
face temperature increases significantly over the range 
of conditions studied. While the dominant, change in 
aerodynamic properties is caused by the change in 




b) 101 km. 



c) 83 km. 

Fig. 3 Number density contours (n/n^,) at various 
altitudes, a = 30°. 
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Fig. 4 Axial force coefficient for the baseline ge- 
ometry, a = 0°. 

density, small variations in the vehicle velocity, at- 
mospheric temperature and composition, and surface 
temperature produce small changes in aerodynamics 
as well. Since the bridging relation only correlates 
the aerodynamic properties to the Knudsen number 
which is essentially just the inverse of the number den- 
sity, changes in the other variables (which are included 
in the DSMC and Navier-Stokes simulations) are not 
taken into account except at the free-molecular and 
continuum endpoints. The sensitivity to trajectory 
variables other than density can be computed at the 
free-molecular condition, and it is estimated that vari- 
ations of up to about. .5% can be attributed to changes 
in velocity and surface temperature along the more 
rarefied portion of the trajectory. These variations are 
sufficient to cause the computed behavior to deviate 
somewhat from sine-squared behavior of the bridging 
relation. While this deviation is not considered signifi- 
cant for the Stardust trajectory, it may be worthwhile 
to investigate alternate bridging methodology for other 
entry trajectories. 

Based on the results in Figure 4, the constants in 
Eqn. 2 were set to a.\ = 3/8 and a.'> = 1/8. These con- 
stants were used to determine the other aerodynamic 
coefficients at non-zero angles of attack with the free- 
molecular and continuum values of these coefficients 
adjusted appropriately. A comparison of the bridg- 
ing relation predictions to DSMC and Navier-St.okes 
results for C'a, C/v and C m)Cg is shown in Figure 5 
for q- = 10°. There is good agreement, between the 
DSMC, Navier-St.okes and bridging relation predic- 
tions of Ca and Cjv, but. the DSMC values tend to be 
higher than both the Navier-St.okes and bridging re- 
lation values for C m)Cg at, the lower Knudsen number. 
At, the Kn=0. 00857, the grid for the DSMC calcula- 
tion was not, resolved to within one mean free path 
near the vehicle surface as it, was for the higher Knud- 
sen numbers. A limited grid sensitivity study showed 
that, refining the grid would reduce the value of C m)Cg 
slightly, but, the resources required to fully resolve the 
Knudsen layer for this flight, condition were not, jus- 


tifiable for a, preliminary design study. Furthermore, 
C m ,cg is very sensitive to small errors in C/v at low 
angles of attack, and the difference in C/v between the 
DSMC and Navier-Stokes results is less than 1% of the 
total force on the vehicle. 





Fig. 5 Comparison of bridging relation predictions 
to DSMC and Navier-Stokes predictions, a = 10°. 

The aerodynamics of the baseline geometry over a, 
wide range of angles of attack is shown in Figure 6. 
The bridging relation is used to provide the complete 
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angle-of-attack results and the DSMC and Navier- 
St.okes results are included for comparison. In general, 
the bridging relation gives a good match to the DSMC 
and Navier-Stokes results for the three aerodynamic 
coefficients shown. Figure 6(c) also demonstrates that 
the vehicle exhibits a static instability in the free- 
molecular regime (Kn=12.8), approaches neutral sta- 
bility at Kn=0.136, and is stable with a trim angle of 
0° near the hypersonic continuum. Although the free- 
molecular result actually has a trim point at slightly 
over 60°, such a large trim angle is undesirable because 
it would require substantially greater thermal protec- 
tion on the afterbody of the vehicle. Furthermore, the 
stability near this trim point is marginal, and the ve- 
hicle could trim in a backward position (o- = 180°) 
where it is statically stable. Even with the addition of 
a reasonable amount of spin stabilization (doubling the 
spin rate from 5 rpm to 10 rpm), six-degree-of- freedom 
simulations have shown that undesirable excursions in 
angle of attack are possible as the vehicle approaches 
peak heating conditions. Therefore, a modification to 
the design was considered that would achieve greater 
aerodynamic stability in the transitional regime. 

Modified Geometry Aerodynamics 

One concept proposed for stabilizing the Stardust 
vehicle in the free-molecular regime is the addition of a 
“skirt’* or disk attached to the afterbody base as shown 
in Figure 2. Free-molecular (“collisionless” DSMC) 
calculations were performed for both the baseline and 
modified geometries over an angle-of-attack range of 0° 
to 180°, and the pitching moment results are shown in 
Figure 7. The addition of the skirt not only provides 
static stability with a trim angle of o- = 0°, but also 
makes the vehicle statically unstable at a = 180° thus 
eliminating the possibility of the vehicle remaining in 
a backward position as it enters the atmosphere. 

However, the presence of the afterbody skirt causes 
undesirable heating and aerodynamic effects in the 
continuum regime. Therefore, a design was considered 
that removes the skirt after the vehicle passes through 
the unstable portion of the transitional regime. A 
bridging relation was constructed that uses the aerody- 
namics of the modified geometry for the free-molecular 
limit and the aerodynamics of the baseline geometry 
(without the skirt) in the continuum limit. In fig- 
ure 8, a comparison of this modified bridging relation 
and the DSMC results for the fully-skirted geometry 
at Kn=0.136 shows good agreement between the two 
and demonstrates static stability with a trim angle of 
zero. 

It should be noted that the use of “Collisionless” 
DSMC to simulate free-molecular flow theoretically 
provides more accurate results than a traditional free- 
molecular analysis based on line-of-sight. techniques^ 
The collisionless DSMC fully accounts for the ther- 
mal velocity spread which can produce higher forces on 


n.QMP. IKr •> si 



a, deg 

a) Axial force, C 


ztn RCMr fk'n-l O ! 



a, deg 

b) Normal force, C'n- 


DRMr, IKn=1 ? RI 



a, deg 

c) Pitching moment, C'm, eg- 
Fig. 6 Aerodynamics of the baseline geometry. 

“hidden” surfaces than the analytical techniques, and 
these forces can become non-negligible at lower Mach 
numbers. Furthermore, the particle-tracing techniques 
used in DSMC account for multiple surface collisions, 
and these multiple collisions may be important for 
concave surfaces such as the skirted geometry. A di- 
rect measure of this effect is obtained by comparing 
the collisionless DSMC axial force coefficient for the 
two geometries at a = 0°. For the unskirted ge- 
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Fig. 7 Pitching moment for the baseline and mod- 
ified geometries in the free-molecular regime. 



Fig. 8 Pitching moment for the bridging relation 
modified to use a combination of baseline and mod- 
ified geometry aerodynamics, Kn=0.136. 

ometry, Ca = 2.063, while for the skirted geometry, 
Ca = 2.112. Based on a simple line-of-sight free- 
molecular analysis for an infinite Mach number, the 
two values would be the same. While this difference 
is not large, it serves to point out the need to account 
for concavities if the most accurate results are desired. 

Flowfield and Surface Properties for Baseline and 
Modified Geometries 

To design an afterbody skirt such that the skirt 
will survive the portion of the trajectory where it is 
needed but will effectively “fail” or be removed due 
to aerothermal loads after passing through the transi- 
tional regime requires detailed knowledge of the forces 
and heat transfer on the skirt itself. Since the skirt 
should survive at least to the point where the unmod- 
ified baseline geometry has neutral stability, it was 
decided to use DSMC to estimate the loads and heat- 
ing on the skirt at Kn=0.136. 

It is instructive to first, examine the differences in 
flowfield structures between the baseline and skirted 


geometries. Figure 9 shows comparisons of the den- 
sity and flow streamlines for both geometries at zero 
degrees angle of attack. The larger wake produced 
by the skirt is expected. However, it is also seen 
that the densities in the concave region adjacent to 
the afterbody are significantly higher with the skirt. 
Although the overall flow is fairly rarefied, the densi- 
ties in the concavity are sufficiently high that there is 
significant recirculation as shown by the streamlines. 
This stagnation-like behavior has a marked effect on 
the afterbody forces and heating as well as on the skirt. 

Figure 10 shows the pressure, friction (shear stress), 
and heat transfer along the body for the baseline and 
skirted geometries as a function of the distance, S, 
along the surface. As expected the skirt has no ef- 
fect on these quantities along the forebody, but pro- 
duces a marked increase in the forces and heat transfer 
along the afterbody. Furthermore, the peak forces 
and heat transfer at the outermost edges of the skirt 
approach magnitudes comparable to those on the fore- 
body. These high forces and heating even at this 
moderately rarefied flow make the task of selecting 
a skirt material and design very difficult. The skirt 
must survive up to and perhaps slightly beyond this 
flight condition, withstand significant heating along 
the outer edges, retain sufficient structural integrity 
to provide the aerodynamic stability, and then be re- 
moved from the vehicle before reaching the hypersonic 
continuum. While the aerodynamic characteristics of 
the skirted geometry are much better than the base- 
line geometry, these design issues coupled with other 
design constraints make this solution impractical at 
the present time. 

Concluding Remarks 

DSMC and Navier-Stokes computations have been 
used to assist in constructing accurate engineering 
methods based on bridging relations to predict the 
aerodynamics of the Stardust Sample Return Capsule 
in the low-density transitional flow regime. Bridging 
relations provide the ability to develop a more com- 
plete database that is needed for six-degree-of-freedom 
trajectory simulations. A simple sine-squared form for 
the bridging relation fits the data computed by the 
more exact, methods reasonably well and meets the cri- 
teria, needed for Stardust, trajectory simulations. Im- 
provements to the overall methodology may be needed 
to provide better representation of the aerodynamics 
in the most, rarefied portion of the transitional flow 
regime where benchmark calculations are required but, 
the present, methodology is adequate for preliminary 
design studies. 

The baseline Stardust, configuration has been shown 
to exhibit, a, static instability in the free-molecular 
regime that, persists well into the transitional flow 
portion of the trajifpt.ory. Eliminating this instability 
poses a, particular challenge because of the high entry 
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